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Abstract. Hydrodynaniic behavior and the associated discussions of viscosity 
at RHIC has inspired a renaissance in modeling viscous hydrodynamics. An ex- 
planation of Israel-Stewart hydrodynamics is presented here, with an emphasis 
(~! on the tangible benefits compared to Navier Stokes. 
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1. Introduction 



Strong collective flow observed at RHIC has validated the use of hydrodynamics to 
model the collision for much of the reaction. However, even with the small inferred 
viscosities, viscous effects are of the order of 10-20 percent. In the last year sev- 
eral significant advances have been made regarding the formalism, implementation 
and understanding of viscous relativistic hydrodynamics for RHIC. A consensus 
(yO has been building that Israel-Stewart treatments are superior to the traditional 

Navier-Stokes approach, both for numerical and physical reasons. In this talk, I 
^ will explain and motivate the Israel-Stewart approach, then review the observable 

manifestations of viscosity in high-energy heavy-ion collisions with an emphasis on 
^ . explaining femtoscopic measurements of source size. 

Viscous effects refer to the alteration of the stress-energy tensor due to the 
system's inability to maintain equilibrium. For equilibrated systems, the spatial 
components of the stress-energy tensor in the fluid's rest frame are: 

T„=P(e,p)%. (1) 

Here, the pressure is a function of the energy density and the conserved charges p. 
For rapidly expanding systems equilibrium cannot be maintained and T*-' differs for 
a variety of physical reasons, [HI H]: 



1. Local kinetic distributions become anisotropic due to anisotropic flows. This 
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is the text book case which leads to the shear viscosity being proportional to 
the mean free path [ |3] . 

2. Interactions and correlations that extend over a finite range contribute to both 
the shear and bulk viscosities [ H] . For interactions of a finite range matter 
at r = will be heated by the collective flow of the matter at distance 

3. Near Tc, mean fields may find it difficult to adjust to rapidly changing equi- 
librium values. This results in peaks to the bulk viscosity near Tc [121 [5] • 

4. Chemical populations, especially for massive particles, fall out of equilibrium 
for rapidly expanding systems, resulting in fugacities. These effects can be 
incorporated either by explicity treating the unequilibrated populations as 
dynamical charges, or as a bulk viscosity. 

5. The system may not have relaxed from initial conditions for T*^ . Most notably, 
at early times the stress-energy tensor might be dominated by longitudinal 
color fields [ El ^ , which are characterized by large transverse components of 
T*^ and a small, perhaps even negative, longitudinal component. 

Historically, viscous effects have been most often incorporated through the 
Navier-Stokes equation, 

= P5,, -CV-v-rj {d^v, + djv, ~ (2/3)V • v) . (2) 

Like the pressure, the bulk and shear viscosities, C and rj, are functions of e and p. 
The coefficients represent fundamental properties of the system that can be related 
to correlations in an equilibrated system through Kubo relations, 

/•oo 

/ dtt{[T^y{r' = 0,t' =0),T^y{r,t)]), (3) 



I dtt{ [5f{r' = 0, t' = 0)(5f (r , t)] ) , 



6f = ^{T,^ + Tyy + T,,)-P-ci{Too~{Too)). 

Lattice Gauge theory provides such correlations of the stress energy tensor, but 
only for imaginary times, < t < if3. The difficulties of extracting to real times is 
discussed in [[Tl]. 

Navier-Stokes treatments can be problematic for two reasons. First, the treat- 
ments can have super-luminar or unstable behavior for short wavelength modes. 
Secondly, the treatments implicitly assume that the anisotropy is determined by the 
instantaneous velocity gradient, rather than the past history or initial conditions. 
This restricts the application to small velocity gradients and to times sufficiently 
long so that memory of the initial state of the stress-energy tensor is lost. Israel 
Stewart approaches [[71[Hl[ni[Tni[ni[ni[T31[Tl] avoid both shortcomings. 
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Israel-Stewart hydrodynamics is predicated on the assumption that the non- 
equihbrium part of Tij will exponentially relax toward Navier-Stokes values, 
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where ai and b are the non-equilibrium parts of the stress energy tensor, 

b = :^{T,. + Tyy + T,,)-P, (5) 
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Here, the exponential scaling functions are not arbitrary, but like the viscosity can 
be expressed as correlations, 

al = J d\{Sf{r' = 0)df{r)), al = j d\{T^y{v' = 0)r,,(r)). (6) 

Furthermore, the relaxation times are related to the fluctuations and to the viscosi- 
ties, 

n = CT/al Ta = iiT/al. (7) 

The constraints in Eq.s (|6| and Eq. ^ follow from the constraints that entropy 
production must always be positive, or in the case of Eq. 0, fr om linear response 
theory [[14l[T5]. Unlike the viscous coeflicients, cr^ and cr^ are equal-time correlations 
and can be calculated from lattice gauge theory without analytic continuation. 
Thus, relaxation times can be determined from lattice with the same confidence as 
viscosity coeflicients. 

Numerical problems with Navier-Stokes hydrodynamics derive from the parabolic 
nature of the equations, which like the diffusion equation has arbitrarily fast modes 
[13 [8]. The odd behavior occurs for high wave numbers and disappears for relax- 
ation times of order, or larger, than a characteristic time r^har = where h is the 
enthalpy density and Cg is the speed of sound [ 14J. If one uses Eq. [7]to estimate 
the relaxation times in terms of the fluctuations of the stress energy tensor, and 
if one then calculates the fluctuation and the enthalpy for independent particles, 
one can calculate both the characteristic time and the relaxation times, which both 
are proportional to 77 for shear waves. The ratio of the two times is independent 
of ?7 and comfortably above the values where singular behavior might occur. Thus, 
if one follows the prescriptions for generating relaxation times consistent with the 
viscosity and fluctuations, the relaxation times should be sufficiently long to avoid 
numerical instabilities. 

Elliptic flow is the most common observable to be associated with viscosity. 
However, viscosity also signiflcantly affect spectra and correlations [[21 [13]. These 
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manifestations can be understood by considering transverse acceleration at early 
times, which is driven by the transverse components of the stress energy tensor, 
Txx and Tyy. Since the velocity gradient is largely longitudinal, the effect of shear 
is to lower the longitudinal component, T^z, while increasing the transverse compo- 
nents. From a variety of perspectives, it appears that an early boost of transverse 
collective flow is needed to explain the i?out/^sidc ratio measured in correlations 
analyses [ 1131 1161 117] . Even though hydrodynamics, even viscous hydrodynam- 
ics, is inapplicable for the earliest times, the equation for transverse acceleration, 
{e + Txx)Dv / Dt = —dTxx, is always valid. In fact, for non-interacting longitudinal 
classical fields, T^x = Tyy — e, three times higher than for an ultra-relativistic gas. 
A non-equilibrium system dominated by longitudinal classical fields is the most 
effective means to jump start transverse fiow. 

In the context of elliptic flow, the behavior of non-interacting flelds contrasts 
with that of non-interacting particles. For non-interacting particles, strong trans- 
verse collective flow also develops, but is accompanied by a canceling anisotropy 
in the stress-energy tensor as measured in the frame of the flowing fluid. I.e., if 
one were looking in the x direction, Txx would be less than Tyy if measured in the 
frame of the fluid. However for classical longitudinal flelds, the stress energy tensor 
remains isotropic in the rest frame, allowing significant elliptic fiow to be gener- 
ated within the first 1 fm/c. Thus, opposite much of the current wisdom, delayed 
thermalization can amplify both elliptic and radial fiow. 

Bulk viscosity is only expected to be significant near Tc [El |5]. This has the 
effect of increasing entropy and slowing the expansion, especially at the surface of 
the fireball. Friction diminishes collective fiow, thus modestly reducing the mean 
transverse momentum. The increase in the entropy is then accommodated through 
increased source volumes as determined by two-particle correlations. 

To better understand these ideas, we have performed calculations based on the 
output of a radially- and boost-symmetric Israel-Stewart hydrodynamic code, which 
then feeds phase space points into a Boltzmann code that models the hadronic stage. 
Using color- glass initial conditions [[201 HI] with an initial stress energy tensor set at 
Txx — Tyy — e/2, and with Israel-Stewart relaxation times of 0.5 fm/c, calculations 
were performed with a simple equation of state that crudely mocks the soft region 
seen in lattice calculations for energy densities ~ 1 GeV/fm'^ and the stiff region for 
higher energy densities where ~ 0.3. Figure [l] compares two calculations, one for 
zero viscosity, and a viscous Israel-Stewart calculation where the initial anisotropy 
was chosen to be Txx = Tyy = e/2, and Tzz — 0. The anisotropics associated with 
the initial conditions and shear reduced the i?out/^sidc ratio by ~ 15%. 

Femtoscopic source radii are affected by a host of model choices and parameters, 
not just the viscosity. Among the features which have a non-neglible impact are: 

• the initial profile - color-glass initial conditions led to modestly lower i?out /-Rside- 

• the equation of state - stiffer equations of state lower i?out /^sidc and decreases 
all three dimensions. 

• adjusting the equation of state in the hadronic state to refiect non-equilibrium 
chemistry - modestly inhibits surface emission and lowers i?out/^sido- 
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Fig. 1. Femtoscopic calculations of i?out/^sido from a hybrid hydrody- 
namic/Boltzmann model for the case with viscosities and an initial anisotropy of 
the stress-energy tensor (circles) match expectations from STAR data (stars) [flS]. 
while the ideal hydrodynamic calculation remains ~15% higher. Similar results 
were seen in [ [T3] 



• relaxing the boost invariant assumption for longitudinal expansion - this low- 
ers i?io„g by 5-10% [ 22 . 

• shear viscosity and early anisotropy of the field-dominated phase - reduces 
i?out/i?side by -20%. 

• bulk viscosity during the mixed phase - increases all three radii. 

The failure of some of the earlier hydrodynamic/Boltzmann models to reproduce 
i?out /-Rside does uot derive from a single cause, but derives from numerous enhance- 
ments and model features, some of which change results by a few percent, while 
others change results by 10-20%. 

It should be emphasized that determining bulk properties from femtoscopic 
analyses alone is naive. Many of the features and parameters from the list above 
also affect spectra and elliptic flow. For instance, the mean pt of protons was 10% 
higher in the viscous calculation, as compared to the non- viscous calculations shown 
in Fig. [l] Given the complex interplay between numerous observables and numer- 
ous features of the modeling, it is unlikely that rigorous quantitative conclusions 
about the bulk properties of matter can be extracted without a global model/data 
analysis. Nonetheless, the tremendous progress in modeling relativistic viscous hy- 
drodynamics during the last year suggests that the model space is sufficiently flexible 
to find at least one point at which data can be reproduced. 
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